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13.  ABSTRACT 


The  behaviour  of  spectra  and  cospectra  of  turbulence  in  the  surface  layer  is 
described  within  the  framework  of  similarity  theory  using  wind  and  temperature 
fluctuation  data  obtained  in  the  1968  AFCRL  Kansas  experiments.  With  approp¬ 
riate  normalization,  the  spectra  and  cospectra  are  each  reduced  to  a  family  of 
curves  which  spread  out  according  to  z/L  at  low  frequencies  but  converge  to  a 
single  universal  curve  in  the  intcrtlal  subrange.  The  paper  compares  these  result.1 
with  data  obtained  by  other  investigators  over  both  land  and  water. 

Spectral  constants  for  velocity  and  temperature  are  determined  and  the 
variability  in  the  recent  estimates  of  the  constants  is  discussed.  .The  high - 
frequenev  behaviour  is  consistent  with  local  isotropy.-  In  the  inertial  subrange, 
where  the  spectra  fall  as  n-5/3.  the  cospectra  fall  faster;  uw  and  w 9  as  n“'/°, 
and  u0,  on  the  average,  as  n"5/2.  The  4/3  ratio  between  the  transverse  and 
longitudinal  spectral  levels  is  observed  at  wavelengths  of  the  order  of  the  height' 
above  ground  under  unstable  conditions  and  at  wavelengths  of  the  order  of  L/ 10 
under  stable  conditions.  This  lower  isotropic  limit  is  shown  to  be  governed  by  f  • 
the  combined  effects  of  shear  and  buoyancy  on  small-scale  eddies.  ,  '  ,  ' 
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Summary 

The  behaviour  of  spectra  and  cospcctra  of  turbulence  in  the  surface  layer  is  described  within  the  frame¬ 
work  of  similarity  theory  using  wind  and  temperature  fluctuation  data  obtained  inthe  1968  AFCRL  Kansas 
experiments:  With  appropriate  normalization,  the  spectra  and  cospectra  are  each  reduced  to  a  family  of  curves 
which  spread  out  according  to  zjL  at  low  frequencies  but  converge  to  a  single  universal  curve  in  the  inertial 
subrange.  The  paper  compares  these  results  with  data  obtained  by  other  investigators  over  both  land  and  water. 

Spectral  constants  for  velocity  and  temperature  are  determined  and  the  variability  in  the  recent  estimates 
of  the  constants  is  discussed.  The  high-frequency  behaviour  is  consistent  with  localusotropy.  In  the  inertial 
subrange,  where  the  spectra  fall  as  n-5/,(  the  cospcctra  fall  faster:  utuand  wO  as  n~in,  and  uO,  on  the  average, 
as  n~m.  The  4/3  ratio,  between  .the'transvcrse.and, longitudinal  spectral  levels  is  observed  at  wavelengths  of 
the  order  of  the  height  above  ground  under  unstable  conditions  and  at  wavelengths  of  the  order  of  L/ 10  under 
stable  conditions.  This  lower  isotropic  limit' is  showri  to  be  governed  by  the  combined  effects  of  shear  and 
buoyancy  on  small-scale  eddies. 


1.  Introduction 


Since  the -early  1950’s,  when  modem  recording  and  computing  techniques  became 
available,  considerable  effort  has  gone  into  the  study  of  atmospheric  turbulence  and  its 
spectral  characteristics.  The  large  amounts  of  data  collected  to  date  by  various  investigators  i 

show  clear  indication  that  spectra  of  wind  velocity  and  temperature  obey  similarity  theory; 
over  a  range  of  frequencies  in  the  surface  layer.  The  1968  AFCRL  (Air  Force  Cambridge 
Research  Laborato-ies)  experiment  in  Kansas  (Haugen,  Kaimal  and  Bradley  1971)  \vas  an 
attempt  to  obtain  a  comprehensive  set  of  data  on  wind  and  temperature  fluctuations  over  a  * 

flat,  uniform  site.  In  this  paper  we  will  use  the  framework  of  similarity  theory  to  describe 
the  spectra  and  cospectra  in  the  surface  layer  over  a  broad  range  of  stability  conditions  and 
to  compare  them  with  results  obtained  byother  investigators. 

Instrumentation  for  the  Kansas  experiment  included  three-axis  sonic  anemometers, 
hot-wire  anemometers  and>  fine  platinum  wire  thermometers  mounted  at  three  levels  , 

(3-66, 11-3  and  22 '6  m)  on  a  32  m  tower.  Surface  stress’ measurements  were  obtained  from 
two  CSIRO  drag  plates  (Bradley  1968)  installed  about  50  and  80  m  to  one  side  of  the  tower. 

Mean  wind  speed  and  temperature  gradients  were  measured  with  standard  cup  anemometers 
and  resistance  thermometers  at  8  levels  on  the  tower  between  2-0  and  32  m. 

Outputs  from  these  sensors  were  sampled;  digitized  and  stored  on  magnetic  tape  by  , 

means  of  a  computer-controlled  data  acquisition. system  (Kaimal,  Haugen  and  Newman  / 

1966).  Analog  signals  from  the  sonic  and  hot-wire  anemometers,  the  fine  pLatinum  wire  ' 

thermometers  and  the  drag  plates  were  sampled  20  times  a  second.  The  hot-wire  anemo¬ 
meter  signals  were  also  differentiated,  low-pass  filtered  (Wyngaard  and  Lumley  1967)  and 
recorded  oh  an  FM  tape  recorder.  Further  details  of  the  site,  instrumentation  and  data 
handling  can  bh  found  in  thd  earlier  paper  by  Haugen  et  al.  ( 1971). 


2.  Data  analysis 

Fifteen  1-hr  runs,  lO.unstable  and  5  stable,  were  selected  for  analysis.  With  data  from 
three  levels  they  comprise  a  iarge:sct  of  spectra  and  cospectra  covering  a  range  of  z/L 
values  from  —2*1  to  +3*3.  The  distribution  of  the  runs  according  to  categories  of  z/L  is 
given  inTable  1.  Mean  profiles,  variances,  correlations  and  other  statistical  parameters  for 
each  of  these  runs  are  presented  in  a  separate  data  report  (Izumi  1971). 
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TABLE  1.  Run  numbers  separated  according  to  zjL  categories.  the  run  number  denotes  the  run; 
THE  NUMBERS  1,  2  AND  3  WITHIN  PARENTHESES  DENOTE  LEVELS  5'66,  11*3  AND  226  m  RESPECTIVELY 


zjL  Range 

Runs  (1  hr  each) 

A 

1 

© 

19(3),  23(3) 

-2-0  to  -1-0 

19(2),  20(3),  23(2) 

-1-0  to  -0-5 

13(3),  14(3),  19(1),  20(2),  23(1),  30(3),  40(3) 

-0*5  to  -0-3 

13(2),  14(2),  20(1),  21(3),  30(2),  40(2),  43(3 

-0*3  to  -0*1 

13(1),’, 14(1),  21(1).  21(2),  30(1),  31(2),  31(3),  40(1),  48(1),  48(2) 

-0*1  to  0 

31(1) 

0  to  +0*1 

37(1) 

+0*1  to  +0*3 

17(1),  24(1),  25(1),  25(2),  37(2) 

+0*3  to  +0*5 

17(2),  17(3),  24(2),  25(3),  37(3) 

+0*5  to  +1*0 

23(1),  24(3) 

+1-0  to  +2*6 

23(2) 

•>  +2-0 

23(3) 

The  spectra  and  cospectra  were  computed using  the  fast-Fourier  technique.  The 
available  bandwidth  was  covered  in  two  stage!).  The  higher  range  (0-005  to  10  Hz)  was 
obtained  by  dividing  each  1-hr  record  into  15  consecutive  blocks  of  4,096  data  points  and 
constructing  a  composite  spectrum  by  averaging  the  16  separate  spectra.  The  composite 
spectrum  was  then  smoothed  by  averaging  spectral  estimates  over  frequency  bands. 

For  the  lower  range  (0-0003  to  0-6  Hz)  a  new  time  series  was  generated  from  'the 
original  series  by  subjecting  it  to  a  16-point  non-overlapping  block  average.  The  spectrum 
computed  from  this  series  has  inherently  more  scatter  than  the  composite  spectrum  but  the 
agreement  between  the  two  in  the  two  decades  where  they  overlap  is  very  good.  The 
higher  range  was,  therefore,  treated  as  our- basic  spectrum  with  estimates  from  the  lower 
range  used  only  to  extend  its  bandwidth  to  0-0003  Hz. 

Friction  velocity  (u#)  values  were  derived  from  the  average  of  two  drag  plate  readings 
as  described  by  Haugen  ct  al  (1971).  Dissipation  rate  (<)  values  were  obtained  from  the 
variance  of  the  differentiated  hot-wire  anemometer  signals  (Wyngaard  and  Cotd  1972). 
Wind  and  temperature  gradients  were  computed  from  the  profile  data  by  differentiating 
curves  fitted  at  5  points,  2  above  and  2  below  the  reference  level  (Businger,  Wyngaard, 
Izumi  and  Bradley  1971). 

The  symbols  used  are  standard  except  where  noted.  V  is  the  mean  horizontal  wind 
vector;  u,  v,  and  w  are  the  fluctuating  velocity  components  along  the  longitudinal  (x), 
lateral  (y)  and  vertical  (z)  directions.  B  denotes  the  mean  potential  temperature  and  8  its 
fluctuating  component  at  any  level.  5  represents  the  average  potential  temperature  for  the 
entire  layer.  Scaling  temperature  T«,  is  defined  as  —u^/u#,  N  is  the  dissipation  rate  for 
0*12  and  k  is  the  von  Kdrmdn’s  constant.  For  reasons  of  symmetry,  this  definition  for  T«, 
is  preferred  over  the  one  with  k  in  the  denominator.  The  dimensionless  terms  used  here  are 


/ 

=■  nzlU 

a  dimensionless  frequency 

zjL 

=  Mg/5)  (Te/u*1) 

a  dimensionless  height 

fa 

=  fez«/u*3 

a  dimensionless  dissipation 
rate  for  turbulent  energy 

4s 

=  kzNlu*TS 

a  dimensionless  dissipation 
rate  for  temperature  variance 

4m 

=  kz(dUldz)lu # 

a  dimensionless  velocity 
gradient 

4ft 

=  kz(d9ldz)IT, 

a  dimensionless  temperature  gradient 
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In  the  discussions  to  follow  we  will  use,  «bois  and  plots  most  commonly  used  in 
atmospheric  work.  We  measure  frequency,  v,t  wavenumber,  spectra;  and  the  conversion 
from  one  to  the  other  is  made  through  the  use  yf  Taylor's  hypothesis.  Strictly  speaking  we 
should  use  different  symbols  to  represent  the  two  types  of  spectra.  Taking  the  u  spectrum 
as  an  example 

j  ~iT‘  =  j  Su(n)dn. 

Jo  Jo 

If  k,  2 wn/U  we  have  i  \ .  > 


SHLw  till] 
u  \uj- 


Su(n), 


,  -  nSit(n). 


<D 


It  is  traditional  in  atmospheric  v »k  to  plot  frequency  spectra,  e.g.  S«(n)  and  n  Su(n), 
not  against  n  but  against  the  nondim  ■  conal  frequency,/.  This  convention  will  be  followed 
in  our  paper.  The  spectral  forms  ->  >olving  the  product  with  n  will  be  referred  to  as 
*  logarithmic  ’  spectra  and  cospecf , 


3.  Evaluation  os  srtial  subrange  spectral  constants 

According  to  Kolmogor  ...  .aw  for  the  inertial  subrange  the  one-dimensional  u 
spectrum  can  be  expressed  as; 

=  •  •  ♦  (2) 

where  is  the  wavenumbcri.ithe  x  direction  (k1  =  2rm/U  by  Taylor’s  hypothesis), 
« is  the  dissipation  rate  and  ir  is  a  universd  constant  estimated  from  various  experiments 
to  be  about  0-51 

The  u  spectra  comp'otea  from  the  Kansas  measurements  are  remarkably  smooth  and 
follow  the  ~r  5/3  power  lew'  with  typical  scatter  no  more  than  a  few- per  cent,  so  in 
conjunction  with  «  values  they  are  useful  for  estimating  oct.  The  values  in  Fig.  .1  arc  based 
on  direct  and  indirect  measurements  (from  hot-wire  data  and  <{>,  curve  respectively)  and 
average ,0-50  ±  0-05  (standard  deviation).  An  earlier  estimate  based  on  a  portion  of  this; 
data  was  0-52  ±  0-05  (Wyngaard  and  Cot6  1971).  Another  analysis  based  on  hot-wire 
spectral  levels  and  <  values, gives  0-53  ±  0-02  (Wyngaard  and  Pao  1972)., 

As  discussed  indeta >! later,  the  inertial  subrange  v  and  tv  spectral  levels  are  higher  than 
the  - u  spectral  levels  by dhe  factor  of  A/3  predated  by  isotropy.  Average  ccj  estimates 
from  v  and  te  spectra  are,  0'48  ±  0-05  and  0-50  ±  0-05  respectively..  Other  recent  experi¬ 
ments  by  Boston  (19Xi)*indicate  a  value  of  0-51  for  while  Paquin  and  Pl  id  (1971) 
report  0-57.  An  unusually  high  value  of  0-69  was  obscr\’ed  by  Gibson,  Stegen  and  Williams 
(1970). 

For  temperature,  Corrsin  (1951)  has  proposed  the  following  inertial  subrange  form: 

F,(k,)  =  fr'-'l’NKr5'3,  .  .  .  (3) 

where  /?,  is  a  constant  analogous  to  «j  in  (i)  and  N  is  the  dissipation  rate  of  9*12.  We  did 
not  measure  N  directly,  but  the  temperature  variance  budget  (Wyngaard  and  Cot6  1971) 
indicates  that  to  a  good  approximation  N  is  equal  to  the  production  rate,  or  —aOidd/dz). 
This  approximation  leads  to  estimates  shown  in  Fig.  1,  which  average  0-82  &  O'08. 
They  agree  well  with  our  earlier  estimate,  based  on  a  portion  of  the  date,  of  0-79  dk  0-10 
(Wyngaard  and  Cot6  1971). 

Unfortunately,  there  are  two  current  conventions  for  the  temperature  variance  dissi¬ 
pation  rate  in.Eq.  (3),  and  this  causes  confusion  when  /Jt  estimates  from  different  experi¬ 
ments  are  compared.  For  symmetry  with  velocity  we  use  the  rate  of  destruction  of  9*12, 

♦Spectra  of  all  three  velocity  component!  are  corrected  for  ipatial  averaging  in  the  ionic  anemomettr  (Kalmal  <1  cl.  1968).  These 
effects  are  Important  only  for  >  ./)  where  I  (the  sonic  path  length)  i!  0'2  m,  but  they  restore  the  -5/3  dope  at  the  high  end  of  the 
spectral  bindwWth. 
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'-2.5  -2,0  -i:5  -l.0‘  -0.5  0  0.5  1.0  1.5  2.0  2.5 

*/L 

Figure  1.  Spectral  contUnti  ec,  and  /?£  plotted  against  tjL.  Direct  <  estimates  are  obtained  from  hot-wire 
measurements  and  indirect  t  estimates  from  the  <j>,  curve  in  Fig.  2. 

calling  it  N.  Other  workers  use  x,  the  rate  of  destruction  of  W1,  and  this  leads  to  a  factor 
of  2  difference  in  &  values. 

Recent  estimates,  as  defined  in  Eq.  (3),  show  considerable  spread.  Panofsky’s 
survey  (1969)  of  Russian  values  shows  a  range  from  0-41  to  0-88.  Gibson  and  Schwarz 
(1963)  measured  0-7  in  the  laboratory  and  Grant,  Hughes,  Vogel  and  Moilliet  (1968)  found 
0-62  in  the  ocean.  More  recently  Paquin  and  Pond  (1971)  obtained  a  value  0-83.  Two  other 
sets  of  estimates  by  Gibson  et  al.  (1970)  and  Boston  (1970),  based  on  direct  measurements 
of  N,  indicate  &  values  of  2-3  and  1-6  respectively.  At  this  point  it  is  not  clear  why>the 
spread  in  /?•  values  is  so  large. 


4.  Spectra  of  velocity  component:', 


In  recent  years  various  efforts  have  been  made  to  bring  together  velocity  spectra  from 
many  sites,  heights  and  thermal  stabilities  (e.g.  Lumley  and  Panofsky  1964;  Bern, an  1965; 
Busch  and  Panofsky  1968)  and  to  define  their  general  behaviour  in  terms  of  similarity 
parameters.  A  number  of  spectral  forms  have  been  suggested,  all  of  which  provide  a  reason¬ 
ably  good  fit  under  near-neutral  conditions;  Our  approach  here  will  be  somewhat  different: 
we  first  collapse  all  spectra  into  universal  curves  in  the  inertial  subrange  and  then  observe 
the  spectral  behaviour  at  lower  frequencies  as  a-function  of  z/L. 

We  note  from  Eqs.  (1)  ana  (2)  that  the  inertial  subrange  logarithmic  tt  spectrum 
normalized  with  u,1  has  the  form 

nS«(n)  .  a,  Jkz(\ll3  /nz\~2li  ... 

ti#2  ~  (2rrk)2'3  W"'-*  \Uj  *  *  '  ’  W 

Using  the  definitions  of  /  and  ^„.this  is 


nSu(n) 

u.2 


(2«fe)ifjJ?*  J 


(5) 


^^CH^GTEIUSTIGSv 


0^5*!aM  aj 0;5,  Eq.  (5)  becomes 

.  .  .  (« 

The  Kansas  wults  haw^veu  the^jliowing  interpolation  formula  for  <f>, ,  based  on  hot-wire 
meisui:eBifcnts  dfi<t  "  ; 

-2<z)L<0  \ 

-  &±4-5\zlL\fl>,  0<z/L<+2/*  ’  •  (7) 


‘^fMSfiajt^ragge:  f  s^^al:;ie^{it'/f*.4  predicted  from  Eqs.  (6)  and  (7)  are 

Fig.  2.  The  agreement  is  good,  with  a  standard 
^Y^o^&wut.a^^t^f^iun^jagr&ment'isTouhd  between  the  0-16  <f>,213  curve* 

t'snd^feiwed'si^c^'develsTppl/iaiiH'iu^ii^^  •  4; ' 


*/L 

•Figure. 2;  .Normaii’ed.iogarithmic  u  spectral  estimates  at/  =  4  compared  with  0-12&l/5  curve  derived  from 

hot-wire  data. 


w? how  itidfide '-fiV*  in  the, normalization  of  u,  v  and  w.spectra  we  remove  the.z/L 
'^P^dcnce  in  their  equstioris.  This  brings  all  spectra,  regardless  of  z/L,  into  coincidence 
m  .the.  inertial  subrange. 

Stirting^th'the  w  s^ctrum,  we  have  the  new  form 


nSu(n ) 

a.2*,2'3 


0-4 


(8) 


A^lot  pfjpg^rithmifctu  spectra  nqhnalized  in  this  manner  is  shown  in  Fig.  3.  The  spectra 
cpitvcrge  to  -ii  ^2/3i)ine=abthevhigh-frequene/  end,  but  at  lower  frequencies  there  is  a 
de|4]r^^^i^#»i;#oording  to  z/L.  The  apparent  smoothness  of  the  spectral 
plots  .arid'tfed^r  .d^af^tipmbetween  different  categories  of  z/L  permit  one  to  draw 

values  of  z/L.  We  thus  obtained  a  family  of 
•catves-represehUng.  Si  -spectre-in  the  range  +2  >  z/L  >  —2.  The  dutves,  seen  in  Fig.  4, 

‘YaKi*  ©f  h  ©ht'&fd  fwoi  ,'Keaus  dtts  (Bwfcaer.et  «!.  1971)  end  u»ed  In  computing  4,.  The  conventional  vsluc  for  *  It  0-4. 
tCWser/hia  Eq.(t)  by  9  fetsr  ©f  //.lisjcjuts  o{  hotropy. 
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.Figure  3.  NorriuiU^logtritJmiic  to  *f»cinti  plotted  »g*ih»t/  *how  low-frequency  end'iepantihg  according 

.to  z/L. 


Figuw'41  Generalized  to  e^AimiSpi  e/L  ysliki  ranging  from  +2*0  to  -  1;0.  StippUng^indicates  absence 

of  <my..\wl!  defined  trend  with  zfL 


in 


the. direction  o^ihcrcaanglx  smaller/ ..as z[L  varies  from  +2*0  to  —0*3.  In  the  range 
*0l3>  z/L  >  *-2*0,  however,  the  spectra  are  not  arranged  according  to  zjL,  but  tend  to 


•ift-.'Sg*..-5i«n#.6.;f;or  z/L>  0„a  and»v, curves  both  display  much  the  same  shape  and 
ibchaviour.as  u>.  "Xhe systematic* progression  with  z/L  observed  on  the  stable  side  breaks 
down  as 
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Figure  5.  Generalized  u.gpectnim  For  z/L.yahtet  ranging  from  -r2’0  to  -2-0.  0+  and  0“ denote  z/L  m  o 
approached  from  iv>itiye  and  negative -aidea  respectively. 


Figure  6.  Generalized  v  spectrum  for  z/L  values  ranging  from  +2-0  to  — 20. 


the  stable  and  unstable  spectra*.  It  appears  as  though  the  spectra  are  excluded  from  this 
region  (indicated  in  these  Figures  by  cross-Iiatch)  and  that  a  sudden  shift  in  the  pre¬ 
dominant  scales  of  motion  occurs  as  z/L  changes  sign.  Longitudinal  vortex  rolls,  plumes 
and  other  convective  circulations  triggered  by  thermal  instability  have  a  profound  effect 
on  the  scale  of  turbulence,  but  the  abruptness  with  which  the  atmosphere  responds  to  a 
change  of  sign  in  the  potential  temperature  gradient  remains  surprising.  The  excluded 
region  in /is  approximately,  an  octave  wide  in  the  u  spectral  plot  and  about  a  decade  wide 
in  the  v  plot.  The  heavy  curve  designated  Ot  in  the  figures  defines  the  neutral  limit  for 
spectra  on  the  stable  side.  The  corresponding  limit,  0",  on  the  unstable  side  is  not  uniquely 
defined  and  does  not  coincide  with  either  of  the  vtwo  envelopes  indicated  by  dashed  curves. 

Examination  of  the  unstable  v  spectrum  reveals  two  distinct  regimes,  one  in  the  range 
/  0‘2  where  it  follows  closely  the  sliape  of  the  neutral  spectrum,  including  the  curvature 

•Tht  time  behaviour  hu  been  reported  by  Dr.  Nitli  E.  Butch  (private  communication)  In  v  ipectra  obtained  at  J’66  m  during 
our  196S  Kaniat  experiment!. 
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near  its  maximum,,  and  another  in  the  range /<  0-2  dominated  by  a  larger,  peak  in  the 
range  0-005  <  /  <  0-02.  These  two  regimes  exist  also  in  the  unstable  u  spectrum,  although 
they  are  not  as  clearly  differentiated  as  in  v.  A  change  in  spectral  slope  at/  c£  0-2  marks  the 
separation  in  the  generalized,  u  plot.  More  discussion  on  the  behaviour  of  the  horizontal 
velocity,  spectra  follows.in  Section  8. 

Comparing  Figs.  4, 5  and  6  with  spectra  reported  by  other  investigators  we  find  good: 
agreement  with  most  measurements  made  in  the  first  100  m.  Fig.  7  shows  examples  of  some 
unstable  spectra  superimposed  on  the  envelopes  defined  by  the  Kansas  data.  The  individual 
data  points  shown  are  spectral  estimates  obtained  over  a  tidal  flat  by  Miyake,  Stewart  and 
Burling  (1970).  Spectra  of  u  and  w  for  Hanford;  Round  Hill,  Cedar  Hill  and  Vancouver  are 
the  composite  curves  published  by  Busch  and  Panofsjcy  (1960).  The.v  spectra  for  Hanford 
and  Round  Hill  are  from  data  reported  by  Elderkin  (1966)  and  Cramer  and  Record  (1969) 
respectively. 


Figure  /.  Envelops*  fey  the.uwtable  Kansas  speeir*  compared  wish  other  recent  sptetws. 


Spectral  estimates  for  the  tidal  fiat  roll  off  more  steeply  at  the  high'frequerjcy.end  than 
other  spectra,  but  this  discrepancj'  is  remPycdwhen.thevestinwt%sre  correctedffpr  Bpafial 
response,  In  the  sonic,  imemoineter.  array;  At.lpw^r  ffequena^  the:  tidal -fiff  :f|«ctra' fit  the 
J^rasxdata  remwkably  y/elllThe  severe  high-frequent  attenuation,  of  die  Hanford.  V 
spectrum  is  probably  related  to  the  time  lag  of  the  seryo,  system- which  hunts;the  azimuth 
wind  direction,  The  Vahcouver;spectra  of  u  and  .derivedffrpmvme^uremenirovcr  the 
•sea  by^Weiler .and  Burlinge(1967), -appear  to  have  slightly  Higher  inertial-subrange  interi* 
sitles  than  five  other  spectra.  .  '  .  ~ 


5.  Spectrum  of  temperature 


form  by  rearranging  term!  in.Eq.  (3).  The  resulting;  expression  is 
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nSi(n)  fa 
T#2  (2rr*)?'3 


(?> 


where  T#  is  the  scaling  temperature  and and  fa  are  the  dimensionless  dissipation  iit<  s 
for  0*/2  and  turb'ilent  energy  respectively.  The  observed  balance  between  the  products  « 
and  dissipation  rates  (i.e.  N'==  —w&d QJdz)  enables  us  to  substitute  fa  for  j>s-  An  empirit  ;il 
relationship  for  fa  has  been  derived  from  the  Kansas  data  by  Businger  et  ah  (1971); 


74(1 -9  z/L)- 
74  + 4-7  z/L, 


i/s 


-2  <z/L 
6  <z/L 


<°  V 

<  w 


(10) 


Since  the  product  d#d«rl/1  Is  a  function  only  of  z/L,  the  expression  in  Eq.  (9)  is  consistent 
with  similarity  theory.  At/  —  4  and  using  fa  =  0-8  and  k  —  0-35,  Eq.  (9)  becomes 


019  fa'  fa~l,s. 


(11) 


The  prediction  in  Eq.  (11)  is  compared  with  observation  in  Fig.  8.  The  agreement  ^  good,, 
with  a  standard  deviation  of  ho  more  than  10  per  cent.  A  roughly  similar  curve  was  i  /btained 
earlier  by  Panofsky  (1969).  Allowing  for  the  difference  in  the  selected/value  and  the  vse  of  k- 
in  the  definition  of-T*,  both  curves  show  intercepts  at  z/L  =  0  which  agree  within  lO  per  cent. 
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Figures.  Normalized  logarithmic  0  spcctril  estimate*  at  J  ■>  4  compared  with  curve  derived 

from  hot-wire  and  mean  profile  data, 


1 


By  including  fafai^in  the  normalization  of  the  logarithmic  temperature  spectrum 
we  collapse  all  inertial  subrange  spectra  into  a  single  curve  (see  Fig.  9).  The  stable 
temperature  spectra,  like  those  of  velocity,  separate  into  distinct  categories  of  z/L  while  all 
the  unstable  spectra  crowd  hrfo  the  relatively  narrow  band  indicated  by  the  stippled  area. 
Even  though  regions  occupied  by  the  stable  and  unstable  spectra  appear  contiguous,  the 
transition  from  one  to  the  other  is  not  smooth.  The  spectrum  shifts  from  0+  to  the  outer 
edge  of  the  stippled  as  z/L  changes  from  positive  to  negative;  with  increasingly 
negative  z/L  thespectrum  moves  back  towards  C4 .  The  trend  is  reversed  when  z/L  moves 
in  the  opposite  direction,  from  negative  to  positive, 
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Figure  9,  Generalized'#  spectrum  for  z/L  values  ranging  from  +2-0  to  -2-0. 


Theprogressive  shift  in  the  0  spectrum  within  the  stippled  area  is  par"  5y  a  consequence 
of  using  T#2  for  norriializatiomfo/T,  has  a  cusp-like  behaviour  near  neutral,  Wyngaard, 
Gotl  and  Izumi  (1971)),  but  the  tendency, of  the  logarithmic  spectral  peak  to  shift  towards 
increasingly  larger/ with  increasing  instability  is  real.,  Lumley  and  Panofsky  (1964)  have 
pointed  out  that  the  location  of  the  fl  spectral  peak  intermediate  between  those  of  u  and  w 
suggests  that  the  fluctuations  in  both  velocity  components  contribute  to  fluctuations  in  0. 
Extending  this  argument  further,  one  could  interpret  the  shift  in  the  logarithmic  0  spectral 
peak  to  larger  values  of/ as  being  the  result  of  a  shift  in  the  relative  influence  of  ti  and  u>  on  0, 
As  the  surface  layerbecomes  increasingly  unstable,  the  influence  of  w  grows  steadily,  while 
that  of  u  declmeo.  This  shift  is  reflected  in  the  trend  of  the  correlation  coefficients  listed  in 
Table  2. 


TABLE  2.  Mean  correlation  coefficient  for  different  categories  of  z/L.  The  time  series  were 

HIGH-PASS  FILTERED  BY  DIFFERENCING  WITH  RESPECT  TO  A  5-MIN  MOVING  AVERAGE 


-=/L 

o 

I 

o 

-0*1  -  0*3 

0*3  -  0*5 

0*5  -  1*0 

1*0 -2*0 

FuO 

—0*56 

—0*41 

-0’iJ 

-0*22 

-0*18 

ruo 

+0*37 

+0*44 

+0*30 

+  6*34 

+0*39 

Tutu 

-0*31 

-0*27 

-0*22 

-0*19 

-0*14 

The  generalized  curves  in  Fig.  9  agree  well  with  most  spectra  obtained  in  recent  years 
(sec  Fig  7).  Conversion  of  spectra  from  other  sources  into  the  format  of  Fig.  7  was  based  on 
the 'observed  <j>N<l>,~lls  behaviour  in  Kansas.  An  exception  is  made  with  the  Bomex  data  for 
which  fa  appears  to  be  five  to\<iix  times  larger  than  for  Kansas  or  other  sites,  but  show 
values  virtually  identical  to  those  from  Kansas.  Here  the  inertial  subrange  fit  was  determined 
on  the  basis  of  /3,  values  reported  by  Paquin  and  Pond  (1971). 

Spectral  estimates  for  the  tidal  flat  (Miyake  el  al  1970)  show  fair  agreement  with  the 
Kansas  spectra  at  /  <  1*0,  but  depart  markedly  at  />  1  *0,  presumably  due  to  noise 
contamination  inherent  in  (their)  sonic  temperature  signals.  Spectra  from  the  other  sources 
follow  the  generalized  inertial  subrange  behaviour  quite  well.  The  composite  curve  for 
Round  Hill  (Panofsky  1969)  does  not  show  the  flatness  at/  or  0-1  seen  in  the  Kansas  spectra, 
but  those  for  Ladner  (McBcan  1970)  and  San  Diego  (Phelps  and  Pond  1971)  show  this 
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clearly.  The  curve  for  Bomex(McBeanl970;  Phelps  and  Pond  1971)  stands  out  in  sharp 
contrast  to  other  spectra;  its  peak  is  shifted  to  a  higher  value  of/ (  ~  0*8)  and  the  spectral 
intensities  on  the  low-frequency  side  are  much  lower.  The  Bomex  curve,  in  fact,  resembles  the 
stable,  spectrum  at  z/L  =  +0-5  in  the  generalized  plot  of  Fig.  9,  -although  conditions 
during  the ''Bomex  runs  are  reported  as  unstable  (-0-33  <  z/L  <  —(Ml).  The  reasons 
for  such  anomolous  behaviour  in  the  Bomex  data  are  still  unclear. 


6.  Gospectra  of  Reynolds  stress  and  heat  flux 

In  this  Section  we  v/ill'examine  the  behaviour  of  the  three  non-zero  cospectra  in  the 
surface  layer:  uw,  tuO  and  u9.  These  are  essentially  the  cbspectra  of  the  shearing  stress; 
vertical  heat  flux  anddiorizontal  heat  flux.  The  cospectrum  is  inherently  more  difficult  to 
measure  than  the  pojver  spectrum  because  the  correlations  between  the  variables  being 
compared  are  sometimes  very  small.  The  cospectra  are  also  particularly  sensitive  tc 
instrumental  errors  which  introduce  '  cross-talk 1  between  tile  two  variables.  For  example, 
instrument  levelling  and  alignment  errors  (Kaimal  and  Haugen  1969)  cause  cross-talk 
between ju^and  w,  while  a  sonic  temperature  signal  has  ^contamination  (Kaimal  1969). 
At  high-wave  numbers,  cospectral  distortion  can  also  arise  from,  spatial  averaging  in  the 
sensor- and  from  separation  distance  between  sensors.  To  avoid  such  distortion  we  have 
limited: Ob'-  /e,  to  less  than1, 1/1,  where  l  (the  sonic  path  length)  is  0-2  m.  The  smallest 
wavelength  resolved  is  then  1-25  m.  In  analysing  o,ur  cospectra  we  adopt  the  same  approach 
jused  with  spectra,  namely,. to  collapse  all  the  curves  into  a  single  universal  curve  at  large/. 
The  functional  relationship  between  the  normalized  cospectral  intensities  at/  =  4  and  z/L. 
are  determined  from  plots^  similar  t>Figs.  2  and  8. 

The  cospectra  of  uw< arid  wO'shovPa  —7/3  power  law  in  the  inertial  subrange.  The 
dimensional  arguments  for  such  a  power  law  and  the  experimental  verification  of  it  from 
the  Kansas  data  are  presented  by  Wyngaard  and  Coti  (1972)  in  a  companion  paper  in  this 
issue.  Other  invesvigators  (Panofsky  and  Mares  1968;  Kukharets  and  Tsvang  1969)  have 
reported  a  — 8/3  cospectral'slope  but  the  more  recent  measurements  of  McBean  (1970) 
appear  to  favour  the  —7/3  slope. 

Cospcctra  of  uO  are  rare  in  the  literature.  The  few  examples  presented  by  Sitaraman 
(1970)  are  limited  in  bandwidth  (/  <  1-0)  and  do  not,  therefore,  provide  any  information 
about  the  inertial  subrange  behaviour  of  uO,  Our  data  show  an  average  slope  of  —5/2  in 
the  range  1  <  /  <  10.  No  dimensional  justification  can  be  made  for  this  power  law,  but  it 
serves  the  present  purpose  of  collapsing  the  cospectra  into  a  single  curve. 


(a)  uw  cospcctrum 

In  the  —7/3  region  the  normalized  logarithmic  uui  cospectrum  should  be  a  function 
only  of  z/L  and/.  Denoting  the  z/L  dependence  by  a  function  G,  we  can  write 

-nCjj^n)^  G(z/L)y-4/Jt 


u. 


where  Cuw(h)  is  the  cospectral  density.  The  function  G(z/L)  determined  empirically  from 
cospectral  estimates  evaluated  at  /  =  4  (sec  Fig.  10)  has  the  form 


1, 

1  +  7-9  z/L, 


-2<z/L^0  ) 
0  ^  z/L  <  +2J ' 


(12) 


To  bring  the  stress  cospcctra  into  coincidence  in  the  —7/3  region  (—4/3  for  +he  logarithmic 
cospectrum)  we  divide  by  G(z/L) 


nCuw(n )  0*56  f_lf, 

u,*G(zlL 3  “  (2»r)</’-/  ' 


(13) 


As  with  spectra  described  in  earlier  Sections,  the  normalization  in  Eq.  (13)  separates  the 
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Figure  10.  Normalized  logarithmic  uto  cospectral  estimates  at  /  *»  4  compared.. with  empirical  formula 

in  Eq.  (12), 


Figure  11.  Generalized  uw  cospectrjm  for  zjL  values  ranging  from  +2-0  to  —  2'0. 


stress  cospcctra  into  clearly  defined  categories  of  z/L  (sec  Fig.  11).  The  position  of  the 
stable  cospectrum  shifts  as  a  function  of  z/L,  while  all  the  unstable  cospectra  crowd  into  a 
narrow  band  indicated  in  Fig.  11  by  the  stippled  area.  The  dashed  curves  approjdmate  the 
upper  and  lower  limits  of  scatter  in  the  unstable  cospectral  estimates.  There  is  some  overlap 
between  the  stable  and  unstable  regions,  0+  being  located  within  the  stippled  area.  This 
overlap  is  consistent  with  the  constraint  imposed  by  our  normalization,  which  requires  the 
integral  of  the  uw  cospectrum  for  z/L  <  0  to  be  unity. 

The  logarithmic  cospectral  peak  is  spread  over  only  two  decades  of  /,  separated 
roughly  as  follow's:  between  0-01  and  0-1  for  z/L  <  0  and  between  0-1  and  1  -0  for  z/L  >  0. 
At  very  low  frequencies  (f  <  0-01),  particularly  at  the  22  6  m  level,  the  cospectral  estimates 
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show  a  tendency  to  reverse  sign  and  become  positive.  Similar  behaviour  has  been* reported 
by  Zubkovsky  and  Koprov  (1969)  and  McBean  (1970).  In  the  Kansas  data  this -sign, 
reversal  is  observed  only  under  unstable  conditions  when  the  low  frequency  cospectral 
estimates  fluctuate  between  large  positive  and  negative  values.  This  erratic  behaviour  in 
the  cospectrum  occurs  precisely  in  the  frequency  range  where  the  logarithmic  spectra  of 
u  and  v  attain  their  maxima,  i.e.,  the  range  where  the  effects  of  surface  layer  convective 
circulations  are  most  strongly  felt.  Some  convective  elements  have  been  fohnd  to  transport 
momentum  upward,  against  the  velocity  gradient  (Kaimal  and  Businger  1970),  others  to 
transport  it  downward  (Haugen  et  al.  1971).  It  is  not  surprising,  therefore,  that  cospectral 
estimates  in  this  range  are  erratic  and  highly  unpredictable. 

Local  free  convection  arguments  (Wyngaard  and  Cot6  1972)  predict  that  G(z/L) 
should  be  constant  under  very  unstable  conditions,  and  as  with  other  similar  predictions 
(Wyngaard  et  al:  1971)  it  is  found  that  this  behaviour  holds  right  up  to  z/L  =  0.  Under 
unstable  conditions,  therefore,  the  inertial  subrange  stress  cospectral  level  is  determined 
solely  by  u,2  and/. 

In  Fig.  14,  the  unstable  stress  cospectra  from  Kansas  are  compared  with  results 
obtained  by  other  investigators.  The  cospectral  estimates  from  Ladner  (McBean  1970) 
agree  well  with  the  Kansas  data  while  those  from  Bomex,  described  in  the  same  reference, 
appear  higher  by  a  factor  of  1*5  at/ >  1*0.  The  composite  cospectrum  for  Round  Hill 
(Panofsky  and  Mares  1968)  falls  off  more  sharply  at  />  1*0,  but  otherwise  shows  fair 
agreement  with  our  results.  The  composite  curves  for  Hanford  and  Vancouver  (also  from 
Panofsky  and  Mares  1968)  show  even  greater  departure  in  the  inertial  subrange,  with 
slopes  and  cospectral  intensities  substantially  different  from  any  observed  in  Kansas. 
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Figure  12.  Norina'frf  logarithmic  utf  cospectral  estimates  at  /  ■«  4  compared  with  empirical  formula 

in  Eq,  (14). 

(6)  to  6  cospectrum 

Applying  the  same  arguments  used  for  the  uw  cospectrum  but  denoting  the  zjL 
dependence  by  a  different  function,  Ht  we  write 
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Figure  13.  Generalized  uw  cospectrum  for  zjL  values  ranging  from  +2-0  to  —2-0. 


H(z/L)  determined  from  the  cospectral  values  evaluated  at/  =  4  (see  Fig.  12)  has  the  form 


/rx  /l.  -2  <  z/L  <  0  \ 

H(z/L)  » +  6<4  zjLi  o  <  z/L  <  +2J  • 


(14) 


The  negative  sign  in  the  proportionality  results  from  the  definition  of  T„;  u#T,  is,  therefore, 
the  negative  vertical  heat  flux.  The  constancy  of  H(z/L)  under  unstable  conditions  is 
consistent  with  local  free  convection  behaviour. 

To  bring  the  w  6  cospectra  into  coincidence  in  the  —4/3  region  we  divide  by  H(z/L). 


nC|ft(n)  1-62  , 4/J 
f  '.T'H(zlL)  (2n)'»J  ' 


(15) 


The  cospectrum  thus  normalized  (see  Fig.  13)  shows  behaviour  very  similar  to  the  stress 
cospectrum,  both  in  regard  to  the  location  of  its  peak  and  the  roll-off  on  the  low-frequency 
side.  But  on  the  high-frequency  side,  the  —4/3  line  is  shifted  at  least  an  octave  higher  on 
the  frequency  scale.  This  extended  contribution  at  the  high-frequency  end.gives  the  wO 
cospectra,  especially  the  unstable  ones,  a  flatter  appearance  in  the  midrange.  As  in  Eq.  (13), 
this  normalization  requires  the  integral  of  the  cospectrum  for  z/L  <  0  to  be  unity,  causing 
the  0+  curve  to  lie  within  the  stippled  area.  The  inertial  subrange  cospectral  level  for  wO 
is  about  three,  times  as  large  as  for  uw.  Smaller  eddies,  therefore,  transport  heat  more 
effectively  than  momentum.  The  ratio  of  the  two  cospectra  in  the  inertial  subrange  can  be 
expressed  as 


CM  H(z/L)T, 

CM  =  G(zlL)u, 


=i  2*9(u$/uu>), 


-2  <  z/L  <  +2. 


(16) 


Since  the  integral  of  the  cospectrum  is  the  covariance,  Eq.  (16)  has  the  corollary  u,  *♦  the 
larger  eddies  (f  <  1*0)  transport  heat  less  effectively  than  momentum. 

The  agreement  between  the  wd  cospectra  obtained  by  McBean  (1970)  at  Ladner  and 
the  Kansas  results  of  Fig.  13  is  extremely  good.  However,  his  Bom  ex  data  show  considerable 
departure;  the  bandwidth  is  narrower  and  the  logarithmic  cospectral  peak  is  shifted  about  a 
decade  higher  than  in  results  obtained  over  land.  A  comparison  of  the  unstable  Tuns  from 
Ladner  and  Bomex  with  the  Kansas  curves  is  given  in  Fig.  14.  The  average  Round  Hill 
cospectrum  (Panofcky  and  Mares  1968)  also  indicated  in  the  Figure,  falls  off  more  rapidly, 
starting  at  least  an  octave  lower  than  the  Kansas  cospectrum. 
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Figure  14.  Envelopes  for  unstable  Kansas  cospectra  compared  with'  other  recent  cospectra. 


(c)  u0  cospcctrum 


The  — 5/2  cospectral  slope  (—3/2  for  the  logarithmic  cospectrum)  indicated  for  u0  is 
not  as  clearly  established  as  arcvthe  —7/3  slopes  for  uw  and  w9.  The  stable  runs  show  a 
slope  of  —5/2  at  all  heights  but  the  unstable  runs  show  slopes  varying  from  —3  (at  5-66  m), 
to  —7/3  (at  22.6  m).  The  theoretical  prediction  of  — 3  by  Wyngaard  and  Cotd  (1972)  is 
approached  only  in  a  few  cases  at  the  lowest  height. 

Using  the  overall  average  of  —5/2  for  the  cospectral  slope  we  write 


nCwi(n) 

u*T# 


oc  K(z/L)/-3'* 


where  K(z/L)  is  a  function  determined  empirically  (see  Fig.  15)  to  be 


1. 

1  +  17-4  z/L, 


-2  <  z/L  <  0  \ 
0  <z/L  <  +2 y 


(17) 


Like  G(z/L)  and  H(z/L),  thisvfunction  h  also  a  constant  under  unstable  conditions. 
As  before,  we  collapse  the  logarithmic  cospectra  at  the  high-frequency  end  by  including 
K(z/L)  in  the  normalization. 


nC««(n)  0-55  „ 

«#T»K(z/L)  “  1 


(18) 


The  normalized  u0  curves  in  Fig.  16  are  spread  over  a  much  wider  range  of  cospectral 
intensities  than  either  uw  or  tu0.  The  stable  cospectra  show  the  same  systematic  progression 
with  z/L  seen  in  earlier  Figures.  The  overlap  between  the  stable  and  unstable  regions  is 
larger  than  for  uw  and  w6,  apparently  the  result  of  using  u#T*  instead  of  u§  for  normali¬ 
zation.  For  z/L  <  0  we  can  use  the  local  free  convection  prediction  of  Wyngaard  et  al 
(1971)  and  write  _ 

u»T*  =  u0/5  <f>m  <j>y,  z/L  <  0.  .  .  (19) 

Although  no  attempt  is  made  in  Fig.  16  to  show  the  trend  within  the  unstable  co- 
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Figure  15.  Normalized  logarithmic  u6  cospectral  estimates  at  /  —  4  compared  with  empirical  formula 

in  Eq.  (17). 
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Fiirure  16.  Generalized  uB  cospectrum  for  z/L  value*  ranging  from  +20  to  —20. 

spectra,  there  is  nevertheless  a  rough  separation  according  to  z/L .  The  top  third  of  the 
stippled  area  is  occupied  predominantly  by  cospectra  in  the  z/L  range  0  to  — 0-3,  the 
middle  third  by  the  z/L  range  — 0-3  to  —TO  and  the  lower  third  by  the  z/L  range  — 1-0  to 
-2-0. 

The  relative  efficiency  of  the  horizontal  and  vertical  heat  transports  in  the  high- 
fvequency  end  of  the  cospcctrum  can  be  estimated  from  Eq.  (15)  and  (18) 

Ctt|(n)  0-34  K(z/L) 

C*,(n)~  (2fl/)l,<  H(z/L)  ’  '  •  •  ^Z0) 

A t,/<=»  1,  the  ratio  varies  from  a  -0-25  (constant  for  z/L  <  0)  to  an  asymptotic  value  of 
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—0-68  on  the  stable  side.  It  is  interesting  to  note  that  the"  ratio  is  considerably  less  than  unity 
in  the  range /  >  1-6,  which  suggests  that  the  smaller  eddies  always  transfer  heat  more 
effectively  in  the  vertical  than  in,  the  horizontal.  The  reverse  is  true  with  the  larger  eddies, 
this  being  most  pronounced  in  stable  conditions  when  u9  cospectral  levels  become  2  —  3 
times  larger  than  the  tad  levels. 


7.  Empirical  formulas  for  neutral  lapse  rats 

The  relationships  derived  in  the  previous  Sections  uniquely  define,  for  all  stabilites, 
the  behaviour  in  the  inertial  subrange.  At  lower  frequencies,  empirical  formulas  can  be 
obtained  for  specific  ranges  of  z/L.  To  enable  comparison  with  other  atmospheric  and 
laboratory  data  we  present  formulae  which  fit  our  neutral  spectra  and  cospectra.  For  w  the 
neutral  spectrum  is  defined  unambiguously  as  the  curve  separating  the  near-neutral  stable 
and  unstable  spectra.  For  all  other  spectra  and  cospectra,  where  the  transition  is  obscured 
by  a  gap  or  an  overlap,  the  limiting  curve  bn  the  stable  6ide,  0+,  is  used  instead. 

It  is  quite  apparent  that  there  is  considerable  similarity  in  the  shapes  of  the  logarithmic 
spectra  and  cospectra.  On  the  high  frequency  side  they  fall  off  according  to  —2/3,  —4/3,  or 
—3/2  depending  on  the  paiameter;  oh  the  low-frequency  side  the  slope  is  very  nearly  +1*0, 
The  empirical  formulae,  therefore,  have  roughly  similar  form. 

At  z/L  =  0,  where  <f>,  =  G  =  H  —  K  ~  1  and  <f>h  =  0-74  we; have 


nS*(n)/«#J  = 

1057/(1  +  33/)5'1  . 

• 

•  (21a) 

nS^nJ/u*1  = 

177/(1'- +  . 

. 

.  (21b) 

nS»(n)/u*J  = 

2//U+ 5-3  (!/)*'»]  . 

. 

•  (21c) 

nS|(n)/T«J  «= 

/53-4//(l  +24/)5's, 
\24-4//(l  + 12-5  /)«>, 

f<  0-15\ 
f>  0-15J  ’ 

.  (21d) 

— nCiac(n)/u*J  == 

14//(1  +  9-6  f)v*  . 

•  • 

.  (21e) 

-nCw,(n)/u#T,  = 

(11  fl(l  +  13-3  /)**», 
\4-4//(l  +  3'8/)2’4, 

f<  i-o) 

5>  i-of 

.  (21  f) 

nCK,(n)/u«T#  = 

4077(1  +  14 f)Vi.  . 

. 

•  (21g) 

The  above  formulae  are  good  approximations  of  the  observed  curves  (see  Fig.  17  for  plots). 
The  only  departure  exceeding  ±  10  per  cent  are  at  the  low-frequency  ends  of  the  vspectrum 
and  the  u0  cospectrum.  Empirical  curves  proposed  by  Busch  and  Panofsky  (1968)  and 
Panofsky  and  Mares  (1968)  are  very  similar  to  those  shown  in  Fig.  17,  except  for  a  small 
displacement  on  the  ordinate.  Although  our  formulae  for  the  cospectra  indicate  asymptotes 
which  are  steeper  than  the  observed  values,  the  equation  fit  the  data  extremely  well  in  the 
range  0-01  <  /  <  4-0. 


8.  Behaviour  at  low  frequencies 

In  earlier  Sections  of  this  paper  frequent  mention  was  made  of  systematic  shirts  in  the 
spectral  and  cospectral  peaks  with  respect  to  changes  on  z/L.. The  most  consistent  behaviour 
in  the  logarithmic  peak  frequency,  fm,  is. observed  under  stable  conditions,  i.e.,  z/L  >  0. 
For  z/L  «*)  only  the  spectra  of  9  and  w  show  a  systematic  trend;  the  others  have  their 
peaks  scattered  over  a  decade  or  more  of/, 

Figs.  18  and  19  shows  plots  of  fm  vs.  z/L  in  the  ranges  where  a  clear  relationship  exists 
between  them.  The  curves  drawn  through  the  data  points  are  fitted  by  eye.  Different 
symbols  are  used  to  identify  the  different  heights  of  observation  and  they  all  apparently  fit 
the  same  curve.  In  fact,  on  the  stable  side,  fm  exhibits  nearly  the  same  trend  with  z/L  for 
all  parameters.  The  following  approximations  can  be  made  for  the  stable  spectral  peaks: 


Figure  18.  Logarithmic  spectral  peak  for  u(  v,  w  and  0  plotted  against  a/L.  Curves  shown  are  fitted  by  eye. 
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(22a) 
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On  the  unstable  side,  fm  for  the  u>  spectrum  continues  to  decrease  with  decreasing  jjt/X, 
but  levels  off  to  a  constant  value  of  about  0-17  when  z/L  exceeds  -1-0.  Spectra  of  ui 
obtained  from  the  Gedar  Hili  tower  (Kaimal  and  Haugen  1967)  show  approximately  the: 
same  value  of  fm  maintained  to  a  height  of  320  m. 

Local  free  convection  predictions,  which  are  successful  for  other  statistics  of  w  and  0 
(Wyngaard  et  al  1971),  should  also  hold  for  their/*  values.  In  lor^d  free  convection  L  loses' 
physical, significance  and  a  becomes  the  only  important  length  scale,  so  that  the  peak 
wavelength  (Am)*,  far  example,  is  proportional  to  z.  It  follow*  that 

[fm)u>  -  2/(Am)»  =  constant, 

as  observed. 

The  peak  frequency  of  the  unstable  temperature  spectrum  shows  the  opposite  trend 
from  w  in  its  approach  to  local  free  convection  state.  Proceeding  from  zero  to  increasingly 
negative  values  of  z/L,  fm  increases  from  O'Ol  to  0-05  in  the  range  0  >  z/L  >  -2-0  and 
levels  off  to  a  constant  value  of  about  0-05  at  z/L  =  -2-0.  The  trend  near  z/L  —  0 
suggests  a  discontinuity  at  neutral  with  a  jump  of  approximately  one  decade  at  the  trans¬ 
ition.  A  similar  plot  presented  by  Panofsky  (1969)  using  data  from  Round  Hill  shows 
greater  scatter,  but  roughly  the  same  behaviour  as  Fig.  18. 

Peak  frequencies  of  u  and  v  also  show  discontinuity  at  z/L  =  0,  the  shift  being  of  the 
order  of  1  to  2  decades.  One  obvious  consequence  of  this  discontinuity  is  the  excluded 
region  if,  the  generalized  spectral -plots  of  Figs.  5  and  6.  The  absence  of  any  dear 
dependence  between  fm  and  z/L  on  the  unstable  side  is  attributed  largely  to  the  influence  of 
mesoscale  features  in  the  atmosphere;  in  fact  there  is  a  tendency  for  n«,  rather  than/m, 
to  be  constant  with  height.  The  constancy  of  tim  with  height  fits  in  with  recent  observations 
(Kaimal  and  Businger  1970)  which  indicate  that  convective  systems  have  vertical  integrity 
so  that  the  time  intervals  between  them  change  little  with  height.  For  the  unstable  u  and 
v  spectra,  nm  falls  invariably  between  0-003  and  0-005  Hz  corresponding  to  time  periods  of 
3  to  5  min,  approximately  the  average  interval  between  the  large  convective, systems  in 
the  surface  layer. 

The  strong  similarity  of  shape  among  the  stable  spectra  suggests  that  a  different 
normalization  might  be  appropriate.  The  normalizations  used  in  the  previous  plots  tended  to 
exaggerate  the  separation  between  peaks  of  the  stable  spectra.  This  vertical  separation  is 
greatly  reduced  by  dividing  the  logarithmic  spectrum  by  the  appropriate  variance*  rather 
than  by  u#*.  A  modified  frequency  scale  ///0,  where  f0  is  the  intercept  of  the  extrapolated 
inertial  subrange  spectrum  with  the  nSa(n)/a*  =  1  line  («  =  u,  v,  w  or  6),  brings  all 
spectra  into  coincidence  on  a  —2/3  line  in  the  inertisl  subrange.  When  all  stable  spectra 
are  plotted  in  this  manner,  they  display  a  universal  shape  (see  Fig.  20)  which  can  be 

*V*rianct  uttd  litre  It  limitttl  to*  bandwidth  0  001  to  10  I U  in  order  to  avoid  th<  influence  oftrtndi  »nd  long-tcrmoKilUtiont. 
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where  <r»;is  the  standard  deviation  of  w.  Since  (/#)»  is,  by  definition,  the/value  when  the 
right-hand  side  of  Eg.  (25)  is  unity 


0  09  (*//»), 


(26) 


where  lu  ==  «r»3/«.  Using  Eq,  (24)  we  have  for  w 

(fmhc  —  zK^m)w  —  0*36  z/l*, 

so  that 


(^m)w  —  3  ltd. 


•  (27) 

•  (28) 


The  wavelength  corresponding  to  the  peak  of  the  logarithmic  w  spectrum  in  stable  con¬ 
ditions  is,  therefore,  three  times  In.  On  the  basis  of  experimental  results  from  other 
turbulent  flows  (Batchelor  1953),  we  can  indeed  expect  lu,  to  be  a  length  scale  characteristic 
of  the  energy  containing  eddies  of  the  w  field. 

Expressions  similar  to  Eq.  (26)  can  also  be  written  for  u  and  v  spectra: 

(So)*  =  ifo'ho  z  0*09  z!U,  ....  (29a) 


where 


(Jo)*  =  (*)V*GV/»-  0*06  z/lti, 
1#  =  °tsl<  and  tic  —  a*} l*. 


.  (29b) 
.  (30) 


Fon  temperature  we  write  from  Eq.  (9) 

nS,(n)  .  fit 

V  (2« )*i»\«*l»«*#*;J  • 


Again  extrapolating  this  behaviour  to  the  nS#(n)/c«2  =  1  line  gives 


(31) 

(32) 


where  1»  ==«M-c>2/Ns/J  is  a  length  scale  for  temperature  analogous  to  the  length  scales 
defined  for  velocity. 

The  stable  logarithmic  cospectra  for  uw  and  wO  also  exhibit  similar  shapes  when  plotted 
in  the  sameimanner  as  the  spectra.  In  Fig.  21,  Jo  is  defined  by  the  intercept  of  the  —4/3 
inertial  subrange  slope  with  the  nCa/j(n)/«d  =  1  line  (a  and  /?  are  the  variables  in  question 
and  a/J  is  the  covariance  in  the  frequency  bandwidth  0*005  to  10  Hz).  A  curve  that  fits 
both  coepectra  is 

nCst(n)  0-mifo) 

%  1  +  l*5(///o)J's  ‘  ’  *  '  K) 

where  3  =?  uw  or  w&. 

The  cospectrumof  uO,  with/0  redefined  by  the  intercept  of  the  —3/2  inertial  subrange 
slope  with  nCut(n)luO  «=»  1,  may  also  be  represented  by  a  single  curve  (see  Fig,  21) 


nCm(n)  0*85(///o) 
ufT~  0=1 1  +  1  •?(///<>)*'*  ‘ 


*  (34) 


Using. cospectral  relationships  expressed  in  Eqs.  (13)  and  (15)  and  (18)  we  can  write 

(foU  =  0*10  (G(zlL)yi* . (35a) 

(foU  -  0*23[H(z/L)]3,<,  ....  (35b) 


(fo)ut  =  Q-U(K(zlL))*l>[~^2l\  . 


(35c) 


The  expressions  for  Jo  in  Eqs.  (26),  (29),  (32)  and  (35)  when  used  in  the  appropriate 
interpolation  formula  give  predictions  of  spectral  and  cospectral  behaviour  in  the  stable 
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Figure  21.  Logarithmic  cospv.ta  of  tiw,  wG  and  uC,  normalized  by  their  covariance,  plotted  *gain»t  modified 
f  sale;  uto  and  utf  iwves  correspond  to  Eq.  (33);  u$  curve  corresponds  to  Eq.  (34). 


surface  layer.  The  unstable  spectra  might  also, be  formulated  in  the  same  manner,  but  the 
shapes  differ  for  each  parameter,  and  more  complex  empirical  formulas  would  be  needed 
(except  in  to)  to  represent  them. 

A  simple  physical  mpdii' explaining  the  variation  of  f0  with  z/L  is  given  in  the 
companion  paper.  This  mmM  (called  local  z-less  stratification)  predicts  that  in  the  very 
stable  limit,  z  becomes  ui&hj  ortant  and  L  is  the  only  significant  length  scale.  A0(  the 
wavelength  corresponding  fojV,  would  therefore  be  asymptotically  proportional  IC’L  and 

fa  —  z/A„  ~  z/L. 

Fig.  5  in  the  companion  p-.p«;  ivows  how  (J0)uw  and  (/0)t«  follow  the  linear  prediction, 
with  Eqs.  (35a)  and  (35b)  }  *ov.diag  the  transition  to  the  predicted  neutral  and  unstable 
behaviour.  Such  linear  trend*  can v  ?  found  in  the  asymptotic  behaviour  of  other/o's  as  well. 

Finally,  we  should  mention  4'  it  the  *  frozen  field  ‘  expression,  Kt  —  2 mjU  (or 
f  —  kx  z/2w)  which  we  have  used  to  convert  from  frequency  to  wavenumber  is  only 
approximate,  particularly  at  low  F  iquencies.  If  the  conversion  were  exact,  a  given  tur¬ 
bulence  field,  whether  measured  from  a  fixed  tower  or  a  rapidly  moving  aircraft,  would 
yield  the  same  spectra  when  potted  against  /.  However,  Lumiey  (1965)  discusses  three 
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reasons  why  ==  2mi/Lf  breaks  down  in' high-intensity  shear  flows:  (i)  the  mean  shear 
causes  aliasing;,  (ii)  eddies  are  continuously  changing  and  are  not  '  frozen’  as  they  are- 
swept  past  the  iristfurnent;.and  (iii)  smaller  eddies  are  imbedded  in  larger  ones,  giving  in 
effect  a  fluctuating  convection  velocity.  He  shows  that  at  high  frequencies  (inertial  range 
and  beyond).6nly  (iii):remainsf.but  this  can  also  be  troublesome  (Wyngaard  and  Pao  1972). 
At  low  frequences,  all-three  become  significant  and  it  is  important  to  note  that  the  effect  of 
the  latter  two  depends  on  turbulence  level  ;(<7u/t/).  The  turbulence  level  of  a  given  field 
appears  lower  to  an  aircraft  observer  than  to  one  who  is  fixed  (the  aircraft  velocity  being 
larger  than  t/),  so  we  expect.aircraft  spectra  to  be  more  nearly  interpretable  as  /^  spectra. 
Spectral  shapes  i  U  the  low / range  would,  therefore,  be  different  for  the  two  observers.  Thus, 
in  the  energy-.containing  region,  Ujn  provides  only  a  rough  approximation  to  the  true 
wavelength. 


9.  Onset  of  local  isotropy 


Although  one  expects  to  find-isot  ropy  at  large /  values,  confirmation  of  this  in  the  first 
20  m  has  been  lacking  (Kaimal,  Borkowski,  Panchev,  Gjosing  and  Hasse  1969).  For 
example,  the-  4/3  ratio  between  the  inertial  subrange  spectral  levels  of  the  transverse  and 
longitudinal  velocity  components  was  not  found  in  earlier  experiments  (Weiler  and  Burling 
1967;  Miyake  et  al.  1970).  Busch  and  Panofsky  (1968)  noted  a  trend  towards  the  4/3  ratio 
in  the‘r.ound  Hill  data;  more  recently  Eidsvik  and  Panofsky  (1970)  have  reported  measure¬ 
ments  over  inhomogeneous  terrain  which  indicate  a  4/3  ratio  for  w  and  u  at  /  >  1  -0.  The 
Kansas  data  show  good  agreement  with  the  4/3  prediction  for  all  but  the  most  stable  cases 
(z/I,  >  1-0).  At  /  =  4,  we  find  the  average  spectral  ratio  to  be  1-33  for  Sw{n)ISu(n)  and 
1-28  for  St{n)ISu{n). 

The  behaviour  of  Sic(n)/Su(n)  as  a  function  of/ is  particularly  revealing  (see  Fig.  22) 
since  it  identifies  the  onset  of  the  isotropic  ratio.  (This  ratio  is  established  at  least  an  octave 
sooner  in  St,(n)/Su(n),  but  its  onset  is  not  as  sharply  defined  under  unstable  conditions.)  To 
avoid  excessive  crowding  of  data  points  on  the  plot,  only  averaged  spectral  ratios  for  the 
different  z/L  categories  (see  Table  1)  are  shown.  Sw(n)/Su(n)  increases  rapidly  in  the 
decade  prior  to  attaining  its  isotropic  value,  moving  systematically  to  higher  values  of  / 
with  increasing  z/L.  Only  the  most  stable  run  clearly  falls  short  of  the  4/3  ratio.  A  unique 
relationship  between  the  u  and  w  spectral  cross-over  frequency,  fe,  and  z/L  is  apparent  in 
the  data.  The  plot  in  Fig.  23  shows  fc  decreasing  monotonically  in  the  range  +2-0  >  z/L 
>  — 2-0,  and  leveling  off  to  a  constant  value  of  nearlyO-2  at  z/L  <  2-0.  The  data  points 
are/«bounded  by  the  peak  frequency  of  the  logarithmic  w  spectrum,  on  the  low 
frequency  side  and  by  (/»)«•,  the  lower  limit  of  the  -5/3  range  in  w,  on  the  high-frequency 
side. 


z/L 

•  <  -2.0 
Q  -2.0  to -0.3 

♦  -0.3  to -0.1 
0  -O.l  to  +0.| 

A  +1.0  to  +0.3 
V  +0.3  to  +0.5 

♦  +0.5  to  +J.0 
o  >  +  1.0 

/--/?//  U 

Figure  22.  Hot  of  the  ratios  of  w  and  u  spectra)  estimates  showing  approach  to  the  4/3  ratio  required  for 

itotropy. 
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Figure  23. 


-2.3  -2.0  -1.5  -1.0  r0.5  0  0.3  1.0  1.3  20 

*/L 


Crossover  frequency,  fe,  for  it  anS'io  spectra  plotted  against  js/L.  Solid  curve  ir  fitted  by  eye"; 
dashed  curves.fedicate  upper  and  lower'  bounds  for  fc. 


At/  >  (Jt)w,  therefore,  ail  three  velocity  components  Have  —5/3  spectra!  slopes  and  the 
spectral  ratios  predicted  by  isotropy.  Furthermore,  at  (ft)w  both  the  wO  and  uw  cospectra 
are  falling  rapidly  (as  the  -7/3  power)  to  zero.  We  can,  therefore,  drop  .the  subscript  w, 
and  treat  ft  as  the  lowe/limit  of  the  locally  isotropic  range.  From  Fig.  23,  fi  is  a  constant 
(-  1-2)  for  z/L  <  Qi  and  approaches  10  z/L  for  z/L  >  0.  The  limiting  streamwise 
wavelength  is,  therefore,  of  the  order  of  z  under  unstable  conditions  and  L/ 10  under  stable 
conditions. 

The  observed  variation  of/i  with  z/L  is  consistent,  as  we  will  now  show,  with  arguments 
put  forth  by  L’thilcy  and  Panofsky  (1964).  Since  anisotropy  results  from  wind  shear  and 
buoyancy, 'they  suggest  that  only  eddies  with  time  scales  small  compared  to  the  time  scales 
of  shear  arici  buoyaucy  will  be  isotropic.  A  suitably  eddy  time  scale  is 

r*=(K  ’£)-»/*  ....  (36) 

where  #  is  the  three-dimensional  spectrum;  it  represents  the  density  of  contributions  to 
turbulent  kinetic  energy;  from  wave  numbers  of  magnitude  n.  We  expect  the  onset  of 
i^trepyto  occur  in  the  *~5li  range  of  E,  so  Eq.  (36)  may  be  expressed  as 

r  =  a ~I/2  t~i/3  K~zli,  ....  (37) 

where#  is  1?  >  three-dimensional  spectral  constant.  As  expected,  this  time  scale  decreases 
with  eddy  d:e. 

There  are  two  buoyancy.timc  scales,  one  associated  with  the  mean  temperature  change 
with  height  and  one  caused  by  temperature  fluctuations.  Including  shear  we  now  have 
three  significant  external  time  scales: 


(dm-1 

\dzj 

(mean  shear) 

.  (38a) 

\g)  I dz j 

(mean  buoyancy)  . 

.  (38b) 

(fluctuating  buoyancy) 

.  (38c) 

Ve  expect  isotropy  for  wave  numbers  where  the  eddy  time  scale  in  Eq.  (37)  is  small 
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compared  to  the  time  scales  in  Eq«(38).  The  criteria  for  isotropy  become 


k 

,,  ,J4  Xn<(dvy‘* 

r->  a-JNjt-J/?,  — 

V  dzf 

.  (39a) 

•  (33b) 

O’ 

*  >  (s|!  y*5’/ 

,  (39c) 

o 

Although  these  criteria  are  derived  for  tl  '<»  hmensionaic  spectra,  they  should-,  be 
applicable  to  one-dimensional  spectra  as  well.  Xn/q.ir  surface-layer  terms,  they  can  be 
expressed  as 

/  >  0-3  (mean.  k,&x)  .  .  .  (40a) 

/>  0-3 

z 

L 

3/4' 

<f>i iV*  (mean  buoyancy) 

.  (40b) 

f>  0-3-j 

j* 

|L 

(fluctuating  buoyant.  ' 

•  (40c) 

Under  very  unstable  conditions  the  buoyancy  criteria  dominate,  and  in  the  i„'Ct>cpnvectios’, 
limit  they  give/ >  0-1.  In  neutral  conditions  only  the  shears  criterion  is  operative,  and  it 
gives  />  0-3.  Under  stable  conditions  all  three  become  proportional  to  z/L;  here,  shear 
is  the  most  restrictive,  giving/  >14  z/L.  Therefor  we  expect  fi  to  approach  a  constant 
under  very  -unstable  conditions,  and  to  yary  as  z/L  on  the  stable  side.  Our  observations, 
Fig.  23,  clearly  support’ this  prediction,  and  suggest  that  we  can  interpret  >  in- our 
inequalitities  as  meaning  greater  by  a  factor  of  5  —  10. 


10.  Conclusions 

The  spectra  and  cospectra  of  velocity  and  temperature  fluctuations  show  systematic 
behaviour  with  z/L  when  expressed  in  appropriate  similarity  co-ordinates.  All  spectra  and 
cospcctra  reduce -to  a  family  of  curves  which  converge  into  single  universal  curves  in  the 
inertial  subrange  but  spread  out  according  to  z/L  at  lower  frequencies.  Other  data  obtained 
in  recent  years  over  land  and  water  fit  the  composite  curves  reasonably  well. 

The  high-frequency  behaviour  is  consistent  with  local  isotropy.  In  the  inertial 
subrange,  where  the  spectra  fall  as  rr5,J,  the  cospectra  fall  faster:  uw  and  wO  as  n"7,J,  and 
u0,  on  the  average  as  rr!,2.  The  4/3  ratio  between  the  transverse  and  longitudinal  velocity 
spectral  levels  is  observed  at  wavelengths  of  the  order  of  z  under  unstable  conditions  and 
about  LI  10  under  stable  conditions. 

The  spectral  constant  aj  agrees  well  with  other  experimental  data,  but  additional 
research  is  needed  to  resolve  the  discrepancy  between  our  (indirect)  f}x  value  and  recent 
directly-measured  values. 
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